Creation, manipulation, and detection of Abelian and non-Abelian anyons in optical lattices 
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Anyons are particle-like excitations of strongly correlated phases of matter with fractional statistics, charac- 
terized by nontrivial changes in the wavefunction, generalizing Bose and Fermi statistics, when two of them 
are interchanged. This can be used to perform quantum computations [ ]. We show how to simulate the cre- 
ation and manipulation of Abelian and non-Abelian anyons in topological lattice models using trapped atoms in 
optical lattices. Our proposal, feasible with present technology, requires an ancilla particle which can undergo 
single particle gates, be moved close to each constituent of the lattice and undergo a simple quantum gate, and 
be detected. 



PACS numbers: 03.67.Lx, 03.65.Vf, 37.10.Jk 

The quest for physical systems where anyons [2] can be 
observed has concentrated so far in effectively 2d materials 
exhibiting topological order [ ]. Abelian anyons, whose inter- 
change generates a nontrivial phase in the wavefunction, exist 
in the Fractional Quantum Hall effect. Non-Abelian anyons, 
whose interchange effects full unitary gates on the wavefunc- 
tion, are expected at certain filling fractions [4] (see recent 
experimental progress in [5]). In spin lattice systems, anyons 
can appear as low-lying excitations of topologically ordered 
ground states (see, e.g., [1,6, 7]). Several implementations 
of lattice models with anyonic excitations have been put for- 
ward [8][9][10][11][12][13][14]. Those involving atoms or 
molecules in optical lattices are especially attractive, given re- 
cent experimental progress [ ]. Specifically, Kitaev's hon- 
eycomb lattice model [ ] can be engineered [11, 12], and 
anyonic interferometry in its Abelian phase can be performed 
with cavity-mediated global string operations [ ] or using 
individual addressing to braid excitations [ ] (but due to the 
perturbative nature of the effective Hamiltonian in this model, 
the visibility of anyonic interferometry is degraded [13, 18].) 

Here we propose a novel scheme to create topologically 
ordered states, generate and braid anyons, and detect their 
statistics for any setup based on particles in optical lattices. 
We use a lattice of particles of species A to build the topo- 
logical code and an ancilla of different species B that can be 
moved independently and brought close to any A particle to 
perform controlled operations on the code [ ] . Preparing the 
ancilla in superposition states, making it interact with appro- 
priate code particles and measuring its state, the following can 
be achieved: creation of a topological state, or a general error- 
correcting code (ECC); creation, braiding, and measurement 
(fusion) of anyons, all operations needed to perform topologi- 
cal quantum computation (TQC) by braiding; and anyonic in- 
terferometry, allowing direct observation of anyonic statistics. 
Note that (i) by using an ancilla with different quantum states 
to perform the manipulation of the anyons, tasks can be car- 
ried that are not possible using classical (e.g., laser) manipula- 
tion of anyons (without single-particle addressability, all pro- 
posed methods lack the power of ours); (ii) there is no need in 



principle of single-particle addressability, specially to perform 
proof-of-principle experiments; (iii) it is based on success- 
fully demonstrated technologies [20] [21] [22]; (iv) it is the first 
realistic protocol for simulating universal TQC in an atomic, 
molecular, optical system system [while engineering the mi- 
croscopic Hamiltonian to build topological protection may be 
some time off (though see [23]), the method herein works in- 
dependently of the existence of the background Hamiltonian.] 

We consider 2d lattices loaded with atoms or molecules, 
e.g., 87 Rb. The ancilla, e.g., 23 Na, can be moved indepen- 
dently using a laser potential not affecting Rb atoms (see 
Fig. 1.) These are now routinely loaded in optical lattices, 
and in the Mott insulator state one can have extended regions 
with one particle per site [22, 24]. Single particles can also be 
loaded in optical potentials and moved without decoherence 
[21]. Our scheme can be extended to layered 3d configura- 
tions [25]. We first consider the toric code Hamiltonian [1], 
with Abelian anyonic excitations only, as a toy model, but our 
scheme is basically model independent; later we apply it to 
the D(53) quantum double model [ ], which has non-Abelian 
anyons and is universal for TQC [26].) 

When the ancilla is brought close to a code atom, they expe- 
rience a 2-qubit unitary Uz = \0) a (0\® I +\l) a (l\®Z between 
their internal levels (X, 7, Z are Pauli operators.) This gate can 
be implemented by cold collisions [27] or any other means 
[28, 29]. Single qubit operations can be applied to the ancilla 
without having to address it, due to the different level struc- 
ture of code atoms. Gates like U x = |0) fl (0| <8>/+ |l) fl (l| <8>X 
can be implemented by applying appropriate gates before and 
after Uz- The internal ancilla state (measurement of Z a ) can 
also be detected with standard techniques without having to 
address it or affect the code. 

We next show how to create Kitaev's toric code [ ] with 
these tools. The code is defined as the ground level of a 
stabilizer Hamiltonian on a square lattice of qubits, realized 
as Rb atoms, at the edges of a square lattice. The Hamilto- 
nian H = — £ V A V — Lp^p is me sum of mutually commut- 
ing stabilizers A v = Yiiev^i anc * B p = YliepZi* wnere v mns 
over all vertices and p over plaquettes and products involve 
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FIG. 1 : 2d optical lattice in the x — y plane loaded with one atom per 
site. Atoms are in different Zeeman levels |0) and 1 1) in their ground 
electronic state, storing quantum information. An ancillary atom of a 
different species in its electronic ground state, with relevant Zeeman 
levels |0) fl and |l) fl , is trapped using laser standing waves along the 
three directions. Lasers in the x — y plane are far detuned from the 
fine structure splitting: The potential controlling horizontal move- 
ment of the ancilla does not depend on its internal state. The laser 
propagating along the z direction is tuned in between the fine-splitted 
excited P levels: The potential controlling the vertical movement de- 
pends on the internal state [27, 28]. 



the qubits surrounding the vertices or plaquettes. We can 
associate the presence or absence of particle-like excitations 
at the plaquettes (magnetic defects) and vertices (electric de- 
fects) with the fulfillment or not of the ground level condi- 
tions A y = +1 and B p = +1. Plaquette and vertex excitations 
are thus characterised by B p = — 1 and A v = — 1 and appear at 
the ends of strings of X and Z operators applied on a ground 
state. These 'particles' turn out to have nontrivial (anyonic) 
exchange statistics due to the anticommutation of the X and Z 
Pauli operators, namely the wavef unction gets multiplied by 
— 1 when a vertex particle winds around a region containing 
a single plaquette excitation. Detection of this phase change 
is possible via interference experiments involving superposi- 
tions of states with and without anyons. Moreover, the de- 
generacy of the code allows to interpret it in terms of a set of 
logical qubits whose Z and X operators are given in terms of 
chains of Z's and X's. 

We work with a rectangular surface with smooth and rough 
boundaries [25], with appropriate 3-body vertex and plaquette 
operators along the boundary providing for a two-dimensional 
code space: One logical qubit is encoded as the eigenvalue 
of a chain of Z's along an edge path connecting the rough 
boundaries. The code space is spanned by +1 coeigenstates 
of the stabilizers. To create a state in the code, we start 
with a well defined state \0)® N (+1 eigenstate of each B v ) 
and measure the A stabilizers sequentially, from left to right 
and top to bottom, using the ancilla. If all outcomes are 
+1, our goal is achieved, since |¥) <* f[ v (l +A V )|0> m [30]. 
If —1 is obtained, we can correct by applying Z b to qubit 
b at the bottom of the vertex using again the ancilla, since 



Z b (l - A v )\0) m = (1 -f A v )|0) w ; Z b is applied to a qubit 
that has not been measured yet (for the last row, Z b can be 
applied to the rightmost qubit.) Once we have measured all 
vertices and thus prepared the state, we can measure all stabi- 
lizers to detect errors and apply error-correcting X's or Z's to 
the corresponding qubits by using the ancilla. We could have 
started at another state, measured all operators in any order, 
and then corrected errors in this way to prepare the desired 
state (this can be used to prepare the target state in models 
beyond Kitaev's.) We now show how plaquette and vertex 
measurements, as well as Z's and Z's, are performed using 
the ancilla. To measure A v , we prepare the ancilla in state 
|+) < 2 OC |0) < 2 + |1) < 2, move it to each qubit in the vertex, and 
apply Ux each time. Then we apply a Hadamard gate to the 
ancilla and measure X a . If the result is ±1, we have applied 
(±|rLev^|+) = (1 ^^v) to the qubits at the vertex, thus 
performing the desired measurement. B v is measured by sub- 
stituting Uz for Ux. To apply X (Z) to a qubit, we prepare 
the ancilla in state \ l) a , approach it to that qubit and apply Ux 
(Uz.) Once the toric code state is prepared, operations within 
the code are performed by applying strings of operators, using 
the ancilla in state \l) a , and applying Ux or Uz sequentially 
on the desired qubits by bringing them close to the ancilla. 
To measure string operators, prepare the ancilla in state |+), 
follow the same sequence, and measure X a at the end. The 
toric code has two kinds of elementary excitations [ ] : pairs 
of anyons in frustrated vertices (electric defects, A v = — 1) and 
in frustrated plaquettes (magnetic defects, B v = — 1), with mu- 
tual Abelian anyonic statistics. They can be created by apply- 
ing Z or X to a given qubit, and can be moved, braided, and 
fused together by applying these operators along a given path 
using the ancilla. Superpositions of states with and without 
vortices, or where they are in different places (see Fig. 2), can 
be created, allowing the observation of fractional statistics: 
The simplest interference experiment is shown in Fig. 3. On 
how to infer anyonic statistics from interference experiments, 
see the Appendix. 

We now outline the preparation and manipulation of anyons 
in a non- Abelian setting universal for quantum computation 
[26]: the lattice quantum double model T>{S^) based on the 
group of permutations of 3 elements S3 [ ] (see a brief discus- 
sion in the Appendix; the full construction is given in [31].) It 
is a lattice model generalizing the toric code, where local de- 
grees of freedom live at the (oriented) edges of the lattice with 
orthonormal bases { \g) } labeled by the six group elements g G 
G. The Hamiltonian, also of the form H = — £ V A V — £ p # p , 
has commuting vertex and plaquette stabilizers imposing con- 
straints on the ground states. Their violations define particle- 
like excitations (anyons) with topological charges (electric, 
magnetic and dyonic), with non- Abelian fusion and braid- 
ing rules. Creation, transport, and fusion of anyons can be 
achieved generalizing the controlled-NOT operations of the 
toric code to controlled left and right group multiplications: 
<' R = |0)b(0| 4 + 1 1)^(1 1 (o L h R )a with a L h \g) = \hg), 
°h\§) ~ Isty' The local degrees of freedom for the D(Ss) 
model are qudits of six dimensions, and their six basis ele- 
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FIG. 2: Top: Creation of a superposition of the vacuum and a pair of 
magnetic defects. On top of the ground state (code atoms shown as 
light circles), an ancilla (dark) is initialised as |+) and brought close 
to a code atom, effecting Ux = |0)(0|®/+|1)(1|®X and creating 
two magnetic defects (big blobs) in the adjacent plaquettes in the 
superposition component where the ancilla is in state |1), while the 
code remains in the ground state in sector |0). Bottom: Anyon trans- 
port. The ancilla interacts via Ux with a code atom between a pla- 
quette satisfying the ground state condition and one of the plaquettes 
containing a magnetic defect in a superposition sector, transferring 
the anyon to the first plaquette. 




FIG. 3: Minimal interferometry experiment for the toric code. Top: 
Ancilla Al, initialised as |1), unconditionally creates a pair El, E2 
of electric defects at neighboring vertices by application of Uz = 
|0)(0| 07+ |1)(1| 0Z. Ancilla A2, initialised as |+), creates a su- 
perposition of the presence and the absence of two magnetic defects 
Ml, M2 at neighboring plaquettes by application of Ux- Bottom: 
Ml is wound around E2 by sequential Ux interaction of ancilla A2 
with the code atoms surrounding E2. Ml and M2 are eventually re- 
annihilated, bringing both sectors to the ground state with a relative 
minus sign, |GS) C0( j e 1 1) Al |— )a2» i- e -> a phase —1 is generated 
in the sector where braiding of defects takes place. In this case, the 
interferometry results can be read from the ancilla lattice by a local 
projective measurement on A2. 



ments can be encoded into ground electronic hyperfine states 
of an alkali atom with enough levels, used as code lattice A. To 
create, transport, and fuse pure electric and magnetic charge 
states (enough to simulate universal TQC [ ]), a 6-state an- 
cilla species B is especially appropriate (see Fig. 4.) As in the 
toric code, vertex operators can be measured using ancilla- 
assisted operations to prepare the ground state (see the Ap- 
pendix.) With one single ancilla, which need not be spatially 
addressed, our algorithm requires 0(nm) steps in an n x m re- 
gion; with an auxiliary lattice with one ancilla per face of the 
code lattice, assuming addressability, it can be parallelized to 
depth 0(n + m) (essentially optimal [ ]; see the Appendix.) 

This scheme is independent of the method used to construct 
the topological state. It can be built by cooling an atomic en- 
semble interacting via an engineered topological Hamiltonian, 
providing in principle topological protection to the code ex- 
cept for anyonic manipulations, which should take the system 
to excited levels in a controlled way (as needed to perform 
TQC as such; on how to simulate relevant Hamiltonians, see 
[11, 12].) But it can also be constructed by the above proce- 
dure using ancillas to impose stabilizer constraints, enough to 
perform proof-of-principle interference experiments, a worthy 
goal by itself. This also allows fault-tolerant quantum com- 
putation with general ECCs (topological codes are excellent 
ECCs: local operators do not mix topological sectors, string 
operators mixing them can be efficiently implemented.) 

The arrangement of ancillas is flexible. One ancilla, in- 
dividually manipulable [ ], suffices in principle. Varying 
degrees of parallelization are possible: the ground state can 




FIG. 4: Level structure for the D(S3) model. The elements of S3 
are encoded into ground electronic hyperfine states of a trapped al- 
kali atom A ( 87 Rb or 23 Na, with 8 levels, or 133 Cs, with 16 levels.) 
A 6-state ancilla B of a different species is used to control opera- 
tions. Ancilla states can be moved independently: bringing \e)s close 
to a code atom A and coupling it to \e)A, a collisional phase gate 
Z e = I a 1b — >b{c\ ® |^)a(^| is obtained; together with simulta- 
neous 1 -particle operations ^allA — ®a v a on the code, it provides 
all controlled operations. Indeed, ^j 1A ^^ a nA yields controlled 
phase gates for any code state with va = \g)A( e \ + k)A(g|> an d trans- 
positions with v A « (\e) + \g})A(e\ + (k) - \g))A(g\ (these can be 
composed to obtain any controlled group multiplication.) 



be constructed with sequential operations on one column of 
the sample (parallel to smooth boundaries) at a time, so as to 
cover the whole sample; in interferometry, commuting opera- 
tions can be done simultaneously; some parallelization can be 
introduced with a coarser optical potential for the ancilla than 
for the code. 

This method shares common problems of optical lattice 
schemes of quantum computation, in particular, spontaneous 
emission. Essentially, only the vertical direction of the ancilla 
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is close to resonance; lifetimes of seconds can be reached by 
tuning the laser between the fine structure levels and can be 
enhanced by restricting manipulation of the vertical dynam- 
ics to the (short) times the ancilla is close to a code atom. 
The ancilla can be repumped after an operation, allowing one 
to repeat the tasks and detect errors. Controlled logic by 
cold collisions [ ] requires cooling the system to the phys- 
ical ground state, possible for both the code (in the 1st Bloch 
band, the ground state of the local potential) and the ancilla. 
Rydberg gates [28, 29] based on dipole-dipole interactions 
eliminate this condition. Ancilla-code interactions must break 
the code in a controlled way, not creating (superpositions of) 
stray anyons spoiling the quantum memory: Theoretical anal- 
ysis [33] and experimental results [34] suggest that excellent 
control and small decoherence rates are achievable. Then the 
implementation benefits from the added protection of topo- 
logical codes (or, in general, ECCs.) With a bias magnetic 
field, arbitrary single qudit unitaries can be realized using fre- 
quency and polarization selectivity of microwave or Raman 
laser pulses. Collisional gates can be realized using trap- 
induced shape resonances [35], or using Raman pulses to map 
code and ancilla ground states {\e)A,B in figure 4) to a vibra- 
tional excited state of each lattice well, evolving by a colli- 
sional phase and mapping back [36]. 

The experimental techniques required by our method are 
(i) independent trapping of two particle species A and B with 
different laser trapping potentials; (ii) diluting the population 
of species B so that each particle is individually addressable; 
(iii) bringing species A to a Mott insulator phase; (iv) initial- 
isation of species A in a product state |O) 0iV ; (v) single par- 
ticle gates on species B; (vi) simultaneous gates on all parti- 
cles of species A; (vii) independent transport of internal states 
of single particles of species B so as to effect cold collisions 
with particles of species A. Additionally, large scale simula- 
tions would require the ability to recool to vibrational states. 
Each one of these techniques has been demonstrated experi- 
mentally; bringing them together will pose an interesting ex- 
perimental challenge. 

A method based on the availability of an ancilla species 
thus allows one to perform all TQC tasks on a given topo- 
logical state, independently of the way in which this state is 
constructed, in optical lattices; e.g., universal TQC based on 
braiding can be performed on top of the ground state of the 
D(S3) model. This scheme is likely to prove the most practi- 
cal and general way to perform TQC in optical lattices. Large 
computations will face a steep scaling (a problem not exclu- 
sive of topological settings), but observing interference phe- 
nomena and applying gates by anyon braiding is feasible with 
today's technology. This method can also be used advanta- 
geously in general ECCs. 
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APPENDIX: SUPPLEMENTARY MATERIAL 

In this Appendix we describe in some detail how to simulate 
non-Abelian anyonic models. Our construction is based on a 
generalization of the toric code which is spin lattice realization 
of the quantum double D{G) of a discrete group G. To real- 
ize this we work on a connected two complex T = { V , jF} 
which is a cellulation of a two dimensional surface described 
by vertex, edge, and face sets. Particles with d = |G| levels 
(qudits) are placed on the edges and physical states live in a 
Hilbert space H = tt^d)®^ where H{d) = {C|j)}^. As 
in the toric code, particles on edges that meet at a vertex v all 
interact via a vertex operator A (v). Similarly, all particles on 
edges that are on the boundary of a face / interact via B(f). 
We pick an orientation for each edge with e = [v/, vjj denoting 
an edge with arrow pointing from vertex vj to v&. The choice 
of edge orientations is not important as long as a consistent 
convention is used. The face orientations are all chosen con- 
sistent with the surface orientation. Whereas in the toric code 
model, the ground states are invariant under local Z2 gauge 
transformations generated by {A(v)}, here we demand gauge 
invariance under transformations: 

T g (v)= f[ L g (ej) [7 R^iej), (1) 

where L g (ej),R g (ej) G U(d), the d dimensional unitary 
group, are the permutation representations of the left and right 
action of multiplication by the group element g G G on the 
system particle located at edge ej. For the particle states 
we make the identification \j) = \gj), where by convention 
|0) = | go) = \e), with e the identity element. The action of 
left and right group multiplication on the basis states is then 
Lh\j) = \hgj),Rh\j) = \gjh). A suitable spin lattice model was 
introduced by Kitaev [1]: 

ffro = -£A(v)-£B(/) (2) 

/ 

where 



{Ue k e3f8~ k ° f{ek) =^0=[vM} 

(3) 

The operator A(v) is a projector onto gauge symmetrized 
states. Similarly, the operator 2fy(v,/), is a projector onto 
states with flux £ at face /. Here the sum is taken over all 
products of group elements gk acting on a counterclockwise 
cycle of edges on the boundary of / such that the product is 
£eG. The function Of {ej) = ±1 according to whether the ori- 
entation of the edge is the same as (opposite to) the face orien- 
tation. Assigning the origin v from which the product is taken 
is important as failure to do so is equivalent to only specifying 
the conjugacy class or magnetic charge. Since the trivial flux 
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is its own conjugacy class we write it B e (v,f) = B(f), which 
is the projection we seek for ground states of Hjo- 

By construction [A(v),A(V)] = [£(/), £(/')] = 
[A (v), /?(/)] = 0. The ground states of //to are then 
manifestly gauge invariant states defined as +1 coeigenstates 
of {A (v), #(/)}. Excited states are described by violations 

of the local constraints A(v),B(/) and are particle-like 

r a i 

corresponding to irreducible representations (irreps) ^ 

where [a] labels the magnetic charge denoting a conjugacy 
class of G, and R{N^) labels the electric charge which 
denotes a unitary irrep R of the normalizer of an element 
in the conjugacy class [a]. Notice that for the group Z2, 
L e = R e = I2 and L gl = R gl = G x , and Hjo is precisely the 
toric code Hamiltonian. 

Henceforth, we fix G = S3 = {gj} 5 j =0 the symmetric group 
of three objects and let r be an n x m square lattice with 
boundary. We place d = 6 qudits on edges labeling each by 
it's location etj-^i = [vij,Vkj] and d = 6 vertex and face an- 
cillas labeled by their locations vj^ and fj^. For a two com- 
plex with boundary the ground state \GS) of Hjo is unique 
[37] and can be constructed by measuring the vertex projec- 
tors in an analogous way to the toric code. We begin with 
all system particles in state \e) = |0) which satisfies the zero 
flux condition (B(f)) = IV/. All vertex and face ancillae 
are prepared in the state |0) where \j) = -4= Lf=o^ 27r 6 |^)* 
We then apply a sequence of operations columnwise begin- 
ning on the left column of vertices {vk,o} and ending on 
the right. For column k we apply the controlled operation 
W(v) = Hhes 3 \h) v (h\® T g (v) between each vertex ancilla Vj^ 
and its edge neighbors. These operations can be done in par- 
allel since [Lh(e),Ry(e)] = 0. Next we measure each vertex 
ancilla Vj^ in the basis {\r}} and given the outcome m(j,k) 
apply the single qudit correction gate Z m ^^ (e/^fc+i), where 
Z(e) = Y, 5 k=0 e 2ink ^\k) e (k\, on the right edge. Finally when 
we reach the rightmost column of vertices {v/ ?m -i} we ap- 
ply a sequence of operations from bottom to top involving 
applying W(v/ jm _i) followed by measurement of the ancilla 
and the correction gate Z m ^' m ~ l \e j^-j^i) on the top edge. 
When we reach vertex v n -\^ m -\ we are done since the oper- 
ator T(v n -\, m -\) is not independent but can be written as a 
product of the others meaning that (A(v n _i j/W _i)) = 1 already. 

Our algorithm has a computational depth of 0(m + n) as 
measured by the number of parallel elementary two qudit 
gates. We find contructions of the operators Y<hes 3 \h)(h\<g>Lh 
or Z heS3 \h) (h\ ®R h in 37 controlled phase gates ^|o)(o|®|o>(o| 
and roughly twice that number of local pairwise basis state 
couplings and suspect this count is optimal. One might won- 
der if a faster ground state preparation procedure is possi- 
ble. The answer is no if the initial state is uncorrected and 
the available set of operations is quasi-local. The reason is 
that the final state has global correlations that are created by 
quasi-local operations. In our algorithm these operations are 
measurements but they could also be adiabatic turn on of the 
summands of Hjo- The time scale to perform the quasi-local 
operations (here the measurement time of A(v)) establishes a 



light cone for the flow of correlations. In [32] it was shown 
by an application of the Lieb-Robinson bound that the min- 
imal time to prepare a topologically ordered state beginning 
in a completely uncorrected state is of the order of the length 
of the correlations. Since the correlation length scales as the 
linear dimension of our lattice, our algorithm is essentially 
optimal. 

The group S3 has three conjugacy classes: the identity 
[e] = {e}, transpositions [t] = {^1^2} and cyclic permuta- 
tions [c] = {c + , c_}, and three irreps: the two one dimensional 
irreps R^ (g) = 1 , the signature representation 

= +l=2?-(c p ), R-(ti) = -l, 

and the two-dimensional irrep 

Ri{e) = h, Ri{tk) = R 2 (c p ) = e'Pf*. 

For D(Ss) there are 8 irreps: the vacuum state II H , pure 

^1 

magnetic charges n^IIyJ, pure electric charges 11^1,11^, 

and dyonic combinations lip j , Ilj^ , 11^ . A complete deriva- 
tion of the fusion rules and braid relations for this model is 
given in [ ]. We focus here pure electric or pure magnetic 
charge states, which in fact are sufficient for universal topo- 
logical quantum computation [26]. 

Pure electric charges are labeled by basis states 
|(/* V) e);(v,-)> where /*, = {§E gSG [/?(g)*W is the 
projection operator onto a subspace belonging to the unirep 
R. Since we are considering pure electric charges, the unireps 
of the normalizer of [e] are equivalent to the unireps of G it- 
self. In the context of the spin lattice model this is interpreted 
as electric charge created by applying the projection operator 
Pjf v onto the system with the group action being local gauge 
transformations 7^(v). Magnetic fluxes are labeled by basis 
states: ,£);(v,/)) and are understood as the result of a 
projection 2fy(v, /). 

Excitations created in the bulk always come in particle 
anti-particle pairs. A generic state of a magnetic charge pair 

a*,*- 1 ]) is 

£e[£] 

with Y,£e[£] \ c £\ 2 ~ 1- The unique vacuum magnetic charge 
pair state invariant under conjugation by fluxes is the state 
with a = l/vT^fW, denoted |0 M ; (v,/), (v'J')). Note the 
vacuum state with neighboring magnetic charge pairs can be 
written 

\0[£}',(vijJij),(vijJij^)) = —= £ Ri{ei-ij-ij)\GS) 

V llfll £e[£\ 

(4) 

We can prepare this by beginning with the ancilla fij in 
the state \0 [£] } fiJ where {|* M ) = Z^O^S" 1 , with Z* = 
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I<£ m ei£} ei2nkm/m K)(£m\ and we have labeled the group el- 
ements in [£] = {/o, . . . and applying the two qudit 
unitary 

F[£](fij) = l|G|-|[fl| ® 1|G| + L ® R i(*i-ijHj) (5) 



Next we measure the face ancilla in the basis {|^). For the 
outcome O^j the target magnetic charge state is created. Oth- 
erwise for outcome , we need a correction step. To do this 
prepare the ancilla fij in the state • and apply the con- 
trolled operation A(vy,/y) where 



A(v,/)=£B g (v,/)®I^(/). 



(6) 



which maps the ancilla / to state \g) f when the flux at / eval- 
uated with base point v is g. Such a controlled operation can 
be decomposed into elementary two qudit controlled rotation 
operators with each edge surrounding / as a control and the 
ancilla as the target. 

In the toric code and in fact for discrete gauge theories for 
all finite Abelian groups [ ], excitations can be propagated 
by applying one local operator which simultaneously annihi- 
lates a charge at one location and create one at a neighboring 
location. This is not true for non- Abelian theories as doing 
so violates a face or vertex constraint. To propagate magnetic 
charges from one face / to an adjacent face f essentially in- 
volves coherently mapping the value of flux at / to the face 
ancilla / using A(v,/) and applying a controlled operation 
on the shared edge of the faces /,/'. After this controlled 
operation the face ancilla / is disentangled from the system 
by mapping the flux at face /' to ancilla f and performing 
a controlled operation between ancillae /,/' and finally re- 
versing the mapping on f . In this protocol we are careful 
to demand only single qudit and nearest neighbor two qu- 
dit interactions. This entire process respects superpositions 
over flux states and can therefore be used to propagate mag- 
netic charges around the lattice [31]. Fusion of a magnetic 
charge pair can be measured by using controlled operations 
to bring the constituent charges in conjugacy class [£] adja- 
cent to one another at faces (/,/') and applying the operator 
A(v,/) = L g (EG^g( v >/) ®L g (f) followed by measurement of 
the ancilla / in the basis |0^ ) /. . . The probability to obtain the 
outcome O^j equals the probability for the pair to fuse into the 
vacuum. 

A generic state of an electric charge pair (R, R*) at vertices 
(v, v') will be represented as an \R\ x \R\ matrix 



|M s ;(v,v')) = -^=£M« v 

\ A „ v 



1*1-1 



L E |(<p, e );(v,-)>K p , e );(v',-)) 

W R-n 



with L|f|=o \M R V \ 2 = 1^1- There is a unique vacuum electric 
charge state which is invariant under conjugation by fluxes: 



1 1^ ; (v, v') ) . The state 1 1^ ; (v/ j , v\ j+i ) ) can be prepared as fol- 
lows. First prepare the vertex ancilla v,-j in state \e) Vij . and 
apply the conditional unitary K(vij,eij ; ij+\) defined by 



*(v,e) = | 



?)e(g\®Rg(v) 



V,* 



L g eG\g)e(g\®R g -i(v) e=[*,v] 



(7) 



Depending on the representation, R, apply a single qubit op- 
eration Wr on the ancilla (v,j): 



Wi 



R 2 



-\e)(e\- 
■ 2\e)(e\ 



\c-)(c- 
-\c-)(c- 



\to)(to\-\h)(ti\-\t 2 )(t2 



(the latter is not unitary but can be constructed using adap- 
tive measurements as detailed in [ ]) and finally apply 
K~ l (vij,eij ; ij+\) to disentangle the ancilla from the sys- 
tem. One can prepare more distant vacuum charge pairs, e.g. 
|l^;(v/j,Vj/j/)) using the same procedure but after each ap- 
plication of K(y,e), swapping the vertex qudit with the next 
relevant vertex qudit on the path (via an intermediary swap 
with the in between edge qudit) then applying the operator 
Ur^Vi/j/) at the end, and finally inverting the steps targeting 
the ancilla. Fusion of electric charge pair \M R ; (v, v 7 )) is mea- 
sured by beginning at vertex v and applying a sequence of 
controlled operations K(v, e) acting on a path of vertices from 
v to v' . At vertex v 7 , the operator L c+ (v f ) is applied and then 
the set of controlled operations along the vertex path is in- 
verted. Finally, the first ancilla v is measured in any unitary 
extension of the basis {\R) = UR^ g es 3 \g)/V6}- The proba- 
bility to measure the outcome state \R.f) is the probability for 
fusion in the vacuum or into the signed irrep. 

We now have all the steps to create and braid and fuse 
anyons. One can verify that all the usual rules for exchang- 
ing and braiding anyons are satisfied [3 ] . In particular it 
is easy to see that braiding charges around each other acts 
trivially since the electric charge pair creation operators are 
diagonal in the logical basis. Furthermore, the gauge trans- 
formations Th(v) can be viewed as creating a magnetic flux 
pair, braiding them around the vertex v and annihi- 
lating. To verify non trivial braiding, one can look for im- 
perfect fusion of electric or magnetic charge pairs into the 
vacuum. Let's describe a simplified interference experiment 
that could be performed in principle on a single face of the 
lattice. We begin by preparing the state of adjacent electric 
charge pairs |1^ 2 ; (vy, v/j+i)) and the vertex ancilla vy in 
where 1^) = (\e) ± \h))/y/2 for some h e S 3 . Next 
we apply the controlled operation W(vjj) which creates the 
state (k)vy|l* 2 ;(vu,v u+ i)> + |A>vy|l* 2 ;(vy,v w+ i)»/^ 
followed by measurement of the ancilla in the basis 1^) with 
outcome m = ± 1 . The outcome distribution satisfies 



P(m= \)-P(m 



-l) = |(i^ ; (v,vO|^W;(v,vO)| 2 

_ |Trfe(/»)]| 2 

1*2 1 2 



which is precisely the fusion probability for R2(h) — > l Rl . 
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It is important to clearly define what we mean by in- 
ferring anyonic statistics from interference measurements. 
For any physical theory, the phase accumulated when 
braiding one particle around another has contributions 
from the statistics fa as well as possible dynamical fa and 
geometric (Berry's phase) contributions fa [39]. Consid- 
ering the interferometry experiment depicted in Fig. 3, 
including the other phases in the process gives the transfor- 
mations \GS)\l) Al \+) A2 -> |£l,£2;(vi,v 2 ))|l) A i|0) A2 + 
|£l,£2)|Ml,M2)|l) Al |l) A2 |£1,£2)|1) A1 |0) A2 + 

^+fe+W|El,E2>|Ml,M2>|l> Al |l>A 2 - |G5)|1>ai(|0>a2 + 
^*+fe+^)|l) A2 )/ v ^. If, e.g., there were a background 
Hamiltonian H = — £/£ v A v — UY* P Bp present, then fa = t4U 
where t is the time to perform the braiding. The statistical 
phase can be extracted by comparing this measured phase to a 
proper reference phase. To do so we can perform a similar set 
of operations but starting in the state | GS) | 1 ) A o | +)a2 such that 
the created electric defects neighboring AO are not enclosed 
by the braiding path of the magnetic defects. Since the 
braiding paths and particle numbers in each component are 
the same in both experiments, the difference in the measured 
phases accumulated on A2 is fa. 

Our scheme is flexible in how the ancillary particles are 
introduced and one could work entirely within a single two 
dimensional lattice with in place face and vertex ancilla. For 
example, the retroreflected square lattice of [ ] provides a 
bipartite square lattice with addressable systems and ancil- 
lary particles. By dynamically changing the trapping fields 
strengths and phases, the face and vertex ancilla can be moved 
along any cardinal direction to interact pairwise with the sys- 
tem particles on the edges. However, to build the ground states 
and manipulate anyons, individual addressability of the an- 
cillary particles using, e.g., the techniques in [40, 41] seems 
necessary. 
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